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SECTION 1 


INTRODUCTION 

In March 1976 Ref. 1 "A User's Guide to Computer Program CIVM-JET 4B 
to Calculate the Transient Structural Responses of Partial and/or Complete 
Structural Rings to Engine-Rotor-Fragment Impact" was published. The general 
capabilities of that two-dimensional (2-D) structural-response finite-element 
computer code are summarized in Table 1. Although large displacements and 
elastic-inelastic behavior of initially-isotropic ductile-metal Bernoulli- 
Euler structures are taken into account, the code is not valid for arbitrarily 
large strains and arbitrarily large rotations for several reasons. First, the 
strain-displacement relation used (see Appendix A of Ref. 1; also called Type B 
in Refs. 2 and 3) is valid for arbitrarily large membrane strains but only for 
small bending strains; the accompanying assumed displacement field employed for 
Bernoulli-Euler behavior is valid also only for small strains since no struc- 
tural thinning is taken into account. In addition, the stress-strain descrip- 
tion used in CIVM-JET 4B is appropriate for small strains but is not valid for 
arbitrarily large strains. Subsequently, the analysis has been extended to 
account properly for large strain behavior of ductile-metal 2-D Bernoulli-Euler 
structures; these developments [2,3] have led to modifications to CIVM-JET 4B 
resulting in the finite-strain computer code now called CIVM-JET 4C. Since 
CIVM-JET 4C is similar in many respects to CIVM-JET 4B + , Ref. 1 can serve also 
as a user s guide for CIVM-JET 4C (thereby avoiding unnecessary duplication of 
elaborate documentation) — together with specific user information included 
in the present report for CIVM-JET 4C. 

Since March 1976, various additions, modifications, and corrections have 
been made to CIVM-JET 4B; most are applicable also to CIVM-JET 4C. In particular, 
in June 1976 some additions (addenda to Ref. 1) were prepared and forwarded to 
the recipients of Ref. 1; these items are contained in Section 2 of the present 
report. Section 3 contains a description of subsequent additions which are 
applicable to both CIVM-JET 4B and CIVM-JET 4C. Next, modifications of the 

+ See Table 1. 
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input and output instructions which are pertinent to CIVM-JET 4C (the finite- 
strain code) are given in Section 4. Additional user instructions for 
CIVM-JET 4C are given in Section 5. 

Finally, it is emphasized that CIVM-JET 4C supersedes and replaces 
CIVM-JET 4B completely since CIVM-JET 4C accommodates finite— strain behavior 
properly for (a) membrane and bending behavior, (b) the material stress-strain 
description, and (c) Bernoulli— Euler deflection behavior including structural 
^hii^riin^ — — whereas CIVM-JET 4B does not account for these features and is 
restricted to "small strains" of undefined size. 
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SECTION 2 


JUNE 1976 ADDITIONS TO CIVM-JET 4B 


In June 1976 the following two additions to the CIVM-JET 4B code of 
Ref. 1 were issued, and are repeated here for completeness; they are applicable 
to CIVM-JET 4C. All cited items have been included in the present versions 
of both codes.* These two additions include; (1) the reactions at structural 
support locations and (2) the translational, rotational, and total kinetic 
energy of each fragment following each impact; these items and the associated 
printed output are discussed in the following under, respectively. Addendum 1 
and Addendum 2. Also added to the printed output immediately prior to the 
printing of the solution data is the following (example) output reminder and 
confirmation of user-specified output options; 


THI POLLOHINR IS THE TIME SOLUTION OF THE FRAGMENT- RING IMPACT. 

OUTPUT HILL BE PRINTED EVERY 40 CYCLES USING OUTPUT OPTION 3. 

REACTION FORCES WILL BE PRINTED AT EACH OUTPUT CYCLE FOR NODES AT MHICH 
B0UN0ARY CONDITIONS AR2 SPECIFIED. D.O.F. THAT APE NOT RESTRAINED AT 
THAT NODE HILL HAVE A REACTION FORCE = 0.0. ALL IMPACTS HILL BE DESIGNATED 
AND ALL THE FRAGMENT ENERGIES HILL BE LISTED AFTER EACH IMPACT. 


Included in Addendum 1 and Addendum 2 are example output data associated with 
these changes for the example problem given in Subsection 6.1 of Ref. 1. 

Addendum 1; Addition of Reaction Force Calculation and Output to CIVM-JET 4B 

For the complete assembled but unrestrained discretized structure 
subject to fragment impact, the equations of motion at time instant t for 
the containment/deflector structure are (pg. 8, Ref. 1) ; 


The program statement sequence numbers shown in this 
CIVM-JET 4B numbers used in Ref. 1. in the CIVM-JET 
statements have resequenced. 


section pertain to the 
4C program, many program 
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( 2 . 1 ) 


where the mass matrix f M*J is a limped (diagonalized) matrix. If the 
structure were supported (that is, certain degrees of freedom i were constrained 
to experience zero displacements, velocities, and accelerations), the governing 
equation for any such ith degree of freedom would read: 


M 




Hence, 


*0 




(2.2a) 


where R*, the reaction force resultant, is the force (or moment) resultant 
applied to the structure by the support; its positive sign is in the direction 
of the defined positive direction of the associated ith pre— constrained degree 
of freedom. 


In order to compute and to print these reaction forces at each printout 
cycle, the following (indicated "added") statements have been added to MAIN 
of the CIVM-JET 4B program: 

In the dimension statement list, added was: 


IMPLICIT REAL*8 (A-B.C-Z) _ ,, . 

DIMENSION CINETF (6) , TBANEN (6) , HOT SB (6) ,HEP(7,3) J aJJfd 

DIME NSICN B5PS (3, 3) , EPI (3) , EPO (3) , HNIN (SI ) 


HAIN0030 
HAZN0032 
MAI NO 040 


Between statements MAIN 5450 and MAIN 5530, the following changes were 

made: 
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686 IP(NBCONDlEQ.O) GO TO 889 
DO 888 I=1,NBCONC 
NXY*NODEB (I) 

IP (NBC (I) . EQ . 1 ) GO TC 686 
IP (NBC (I) .BQ. 2) GO TO e87 
IP (NBC (I) .50.3) GO TO 885 
886 REP (1 , 1) = FLVA ( NXY * 4-3) "1 

PEP (1, 2) = 0.!) zJdei 

REP(I,3) * FLVA (HXX*4- 1) 

PLVA (NXI»4-3) = 0.0 
plva nxi*4-i) =0.0 

GO TO 888 

897 HEP (1,1) = FLVA (NXY*4-3) 

HEP (I> 2) = PLVA (NXY*4-2) 

REP ( I, 3) = PLVA (N X 1*4- 1) 

PLVA (NXY*4- 3) - O.C 
PLVA (N X 1*4-2) =0.0 
PLVA (N XY*4- 1) =0 . 0 
GC TO 888 
885 REP (I, 1) = PLVA (HXY*4-3) | ,, , 

REF(I,2) = PLVA (NXY*4-2) 

REF (1 , 3) = 0.0 -J 

PLVA (NXY*4-3) =0.0 
PLVA (M X Y* 4-2) =0.0 
888 CONTINOE 



MAI N5400 

HAIN5410 

HAIN5420 

HAIN5430 

HAIN5440 

H AIN5450 

HAIN5451 

B AIN5452 

BAINS453 

HAIN5460 

MAIN5470 

HAIN548C 

HAIN5481 

MAIN5482 

HAIMS 483 

HAINS490 

MAIMS 500 

MAIN5510 

HAIK5520 

MAIH5521 

BAIN5522 

BAIN5523 

BAIN5530 

HAIN5540 

MAI NS 550 


Between FORTRAN cards MAIN 7780 and MAIN 7790, the following FORTRAN 
statements have been added: 


7182 FORMAT ( • • , 4X , 1 3, 16 X, CIS. 6, 7X, 13 , 6 X .15, 5X . D 1 5 . 6) 
IF (NBCOND.EQ.O) GOTO 1254 
WRITE ( MWFITE, 1256) 

1256 FORMAT ('C',6X, 'REACTIONS AT NODE* , 15X. • RV (LBS) • . 
allX, 'PM (IN-LES) •) 

DO 1252 1 = 1 , N BCOND 
NX Y= NCCEE(I) 

WRITE (HWRITE, 1253) NXY, (REP(I,J) , J=1,3) 

1253 FORMAT (• ',1 8X , 14 , Sx , 3C2 0 .6) 

1252 CONTINUE 

1254 CONTINUE 
WRITE(MWPITE, 11100) 


13X, ' RW (LBS)' 


hi 

$ 

<o 


HAIU7780 

HAIR7781 

HAIK7782 

MAIH7783 

HAIN7784 

HAIH7785 

HAIH7786 

HAIN7787 

BAIM7788 

HAIN7789 

BAIM778A 

HAIM7790 


These modifications will print out the reaction force resultants RV 
(tangential, lbs) , RW (normal, lbs) , and RM (moment, in- lbs) for each constrained 
node at each user-designated printout cycle. The physical positive orientation 
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of each reaction force (applied to the structure) is the same as the orienta- 
tion of the corresponding local degree of freedom depicted in Fig. 3 of Ref. 1 
(page 224) . Note that the reaction force is set equal to 0.0 for each uncon- 
strained degree of freedom at a "constrained" node; see FORTRAN statements 
MAIN 5452 and MAIN 5523 above. 

These modifications result in the following output for the Subsection 6.1 
example of Ref. 1; the output occurs for cycle 1100 at 132.2 microseconds after 
initial' impact (see pages 171 and 189 of Ref. 1) : 


REACTIONS AT NODE 

RV(LES) 

RW (LBS) 

RH (IN-LBS) 

1 

0* 162271 D+C<J 

-0.595815D+01 

0.0 

13 

0. 157820C+03 

-0. 1491600*04 

0.292786D+04 


******************************************************************************* 

Addendum 2; Printout of Kinetic Energies of Every Fragment After Each Impact 

For every fragment after each impact, the kinetic energies in inch-pounds 
translational (TE) , rotational (RE) , and sum or total (TOE) are printed out. 
This has been accomplished by adding the indicated statements between FORTRAN 
cards MAIN 5630 and MAIN 5640. The correct dimension statement MAIN 0032 of 
Addendum 1 has been included; variables CINETF(6), TRANEN(6), and R0TEN(6) 
from that dimension statement apply to this addendum. 


If (IMCO.EQ.O) GOTO 527 
IHC0 = 0 
DC 6301 1=1, HF 

TRANEN (I)= FHASS(I)/2.0* (UD0T (I) **2*SC0T (I)* *2) 
R0TEN ( I) = FMOI(I)/2.0 * AE0T(I)**2 
CINETF (I) = TRANEN (I) *F0TEN (I) 


J = I-2 
GOTO 6307 

6303 J= J-2 

6307 IF (J) 6304,6305,6303 

6304 VRITE(RNRITE,6302) I, TRANEN (I) ,R0TEN(I) ,CINETF(I) 
6302 FORMAT (• FBAG , ,I2, , J IE*’, D12.4,' RE a *, 012.4,' 


TOE*' ,D12..4) 


GOTO 6301 , 

6305 WRITE ((WRITE, 6306) I, TRANER (I) , R0TER (I) , CIHETF (I) 

6306 FORMAT (• ♦ ' , 62X, • FRAG* ,12,' : TE=',Dli.«,' RE*',D12.4,« 


0012. 4) 

6301 COHTINUE 

00 526 0*1 , NS 


TOE*' , 


A) 

■v 

« 


HAINS620 
HAIN5630 
MAIN5631 
MAIN5632 
MAIN5633 
MAIN5634 
MAIN 56 35 
HAIN5636 
NAIN5637 
HAIN5638 
HAIN5639 
HAIN563A 
M AIN563B 
MAIN563C 
MAI N563D 
MAIM563E 
HAIN563F 
MAIN5640 
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An example of this added output is shown as follows for the example 
problem in Subsection 6.1 of Ref. 1 (see pages 171 and 186): 


IMPACT NO. 1 TINE 0.967796D-03 DUFING CYCLE 968 ELEfl 9 

FRAG 1: TF» 0 . 14650+04 RE= 0.0 TOE* C.1465D+04 
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SECTION 3 


SUBSEQUENT ADDITIONS APPLICABLE TO BOTH 
CIVM-JET 4B AND CIVM-JET 4C 

First, note the following clarifications pertinent to Ref. 1 (with the 
page number of Ref. 1 cited) : 

TPRIM(I) must be a multiple of the At actually used. 

The (circumferential*) "strain" (such as SI or SO) is denoted by 
£ (or y ) and is identified by Eqs. A. 12 through A. 14. 

In Subroutine IMPACT, statement IMPA0720 has been changed to include 
a print message before CALL EXIT to remind the user that the number 
of impacts occurring during a single At time step has exceeded the 
specified 50 "allowable impacts". 

Note that statements ROOT0510 and ROOT0570 in Subroutine R00T4 may 

need to be changed if your computer does not accommodate numbers up 
+50 

to 10 ; the user should change accordingly. 

The quantity e in Eqs. A. 50 and A. 51 is the coefficient of restitu- 
tion: perfectly-elastic impact e = 1; perfectly inelastic impact 

e = 0; and intermediate 0 < e < 1. 

The following "subsequent" additions (since June 1976) are identified 
for convenience and consistency as Addendum 3 through Addendum 7; pages cited 
refer to Ref. 1. All of these items have been included in the present versions 
of both codes + . 

* 

Since this report deals with both initially- curved and initially-straight 2-D 
structures (rings and beams, respectively) , the term "circumferential strain" 
is used throughout this report to denote (a) circumferential strain for rings 
and (b) spanwise strain for beams. 

+ Footnote on page 3 applies. 



118 


153 


257 
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Addendum 3: User Specified Impact-Affected Length L 

eff 

The impact-affected length L* ff (named EFLN(L) in the program) in inches 
may be specified by the user, or the value AtlE/p^ 1 ^ 2 described on page 59 
may be chosen as a default option. Accordingly, the following changes have 
been made, and the new input card requirements are: 

Card 5A on page 32 should now read: 

NSF(L) ,B(L) ,DENS(L) ,DS(L) ,P(L) ,EFLN(L) with format I5,5D15.6 


Card 6 on page 36 should now read: 

DELTAT,DS (1) /P(l) /NTOVR/EFLNf 1) with format 3D15 .6 ,15, D15 . 6 


If EFLN (I) for 1-1/ 2, ... NBR+1 is input as zero or as a negative number , 
the program will use the default value At [e/P q ] 1/ ' 2 . 

In the MAIN program, replace statements MAIN1210 and MAIN1220 with 
the following three cards: 


BEAD (BREAD, 2) DELTAT, DS (1) , P (1) , NTOVR,EFLN ( 1 ) , (EPS (1 ,L) , SIG (1,L) , 

1 L — I 4 ) 

2 FORMAT (3D15. 6,15, D15. 6/(4015.6) ) 


MAIN 1210 
MAIN1215 
MAIN 1220 


Also,i replace card MAIN4690 by the following three cards: 


IF (EFLN (I) .GT. 0.) GO TO 5551 
EPLN(I)=(YOUNG (I) /DENS (I) ) **0.5*DELTA1 
5551 CONTINUE 


MAIN4635 

MAIN4690 

MAIN4695 


In Subroutine Bran, replace BRAN0380 and BRAN0390 with the following two 
cards : 


B£AD(MSEAD,55C0) 
5500 Por..1AT(I5,5D15,6) 


NS Ft (IB) ,B(IB) 


, DENS (IB) , DS (IB) ,? (IB) , EFLN (IB) 


BRAN 03 80 
BRAN0390 


See also Subsections 2.2 and 2.5.2 of Ref. 2. 
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Addendum 4: Addition of Clamped-Sliding Nodal-Displacement Condition 


The following third type of prescribed-displacement nodal condition 
should be added to page 4 of Ref. 1: 

(3) Clamped-Sliding (w = = 0) 

Card 14B on pages 41 and 42 

Use NBC (I) = 4 for clamped-sliding condition; will set w = 0 

and i(i= 0 at NODEB(I) . 

In the MAIN program after MAIN5450, add the following card: 

IF (NBC (I) . EQ. 4) GO TO 8855 


In MAIN after MAIN5540, insert these six cards: 


GO TO 888 
8855 REF (I, 1) =0. 

REF (1,2) =FLVA (NXY*4-2) 
REF (1,3) = FLVA (NXY*4-1) 
FLVA (NXY*4-2)=0. 

FLVA (NXY*4-1)=0. 


MAIN554 1 
MAIN5542 
MAIN5543 
MAIN5544 
J1AIN5545 
3AIN5546 


In Subroutine IDENT after IDNT0910, insert: 


' IF (NBC (I) ,EQ, 4) WRITE (8 WE ITS, 24) NODEB (I ) IDNT091S 


In Subroutine IDENT after IDNT0950, insert: 

24 FORMAT {' CLAMPED SLIDING DISPLACEMENT CONDITION AT NODE =',I5) IDNT0955 
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In Subroutine QREM replace QREM2710, QREM2720, and QREM2730 with these 
three cards: 


IF (NBC (I) • NE. «) CALL EEC (JT433, SPRIN, NI,ICCL) 

IF (NBC (I) .NE. 3) CALL EEC (JT4H1 , SPRIN, NI, ICOL) 
IF (NBC (I) . NE. 1) CALL EEC (JT4H2, SPHIN,NI, ICOL) 


QREf'ZTiP 

QREM2720 

QREM2730 


In Subroutine TSTEP, replace TSTP0480 through TSTP0530 with: 


IF (NBC (1) . HE. 4) 
IP (NBC (I). NE. 4) 
IP (NBC (I) . NE . 3) 
IF (NBC (I) . NE. 1) 
IF (NBC (I) . NE. 3) 
IP (N DC ( I) .NE. 1) 


CALL EEC (JT4a3,STIFK,NI, ICOL) 

TRIAL(JT4M3) - 0.0 

CALL EEC (J”4a1,STIFK,NI,ICOL) 

CALL ERC(JT4H2,3TIFK,NI,IC0L) 

TRIAL (JT4M1) -0.0 

TRIAL (JT4M2) =0.0 


T5TP048C 
TSTP0490 
TSTP0500 
TSTP05 10 
TSIP0520 
TSTP0530 


Addendum 5: Definition of Variable NS in Subroutine IMPCTE 


The following two cards should be added to subroutine IMPCTE after 
statement IMPT0200 to correct an omission: 


ns=IK 

- IF (ICP ,GT . 0) NS=IK*1 


IMPT020 1 
IMPT020? 


Addendum 6: Additional Strain and Relative Elongation Output 

Heretofore the circumferential strain and relative elongation have been 
computed and provided as output as, respectively, SO and EO at the outer surface, 
and as SI and El at the inner surface of the structure at locations designated 

as "additional strain points" — see page 50 of Ref. 1. Added now is the 

. o2 °2 

circumferential strain (called SM) and relative elongation (called EM) 
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at the midsurface ( C° = 0) at each "additional strain point". Accordingly, 
the following changes are needed (have been made in the CIVM-JET 4C source 
code and tape) : 

The cards sequenced MAIN7080 and MAIN7090 were changed to read: 

8707 FORMAT ( 'STRAIN AT ADDITIONAL POINTS' ,9X, 'SI' ,14X, 'SO' ,14X, 'SM' , 

@14X, 'El' ,14X, 'EO' , 14X, 'EM') 

After the card sequenced MAIN7550, insert: 

EM=DSQRT(1.0+2. *FARE) —1.0 MAIN7551 

The card sequenced MAIN7560 was changed to read 

WRITE (MWRITE, 8781) IS, EPASI,EPASO,FARE,EI,EO, EM MAIN7560 

The card sequenced MAIN7570 was changed to read 
8781 FORMAT ( ' ’ ,10X,I3,16X,6D16.7) MAIN7570 

Addendum 7: Corrections to Impact Calculations 

In Subroutine IMPCTE after statement IMPT0180, insert: 

COMBOS /COO/ IB COO IMPT0185 

In subroutine IMPCTE after statement IMPT3300, insert: 

IBCOO=IBCOO-1 IKPT3305 

In Subroutine IMPACT modify statement IMPA0060 to read: 

2NTSD (6) ,NEF (6), EL (50) ,EFLN(6) ,EPCC[51,6) IBPA0060 

In subroutine IMPACT after statement IMP AO 300, insert: 

DDBLBs 0.1D-16 IHPA0305 


MAIN7080 

MAIN7090 
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In Subroutine IMPACT after statement IMPA0340, insert: 


dJJ=0 


IMPA0345 


In Subroutine I IMPACT after statement IMPA0760, insert: 


JF1 = JF 
PALI =P AL 
I BIG1 = I BIG 


IMPA0761 

IMPA0762 

IMPA0763 


In Subroutine IMPACT after statement IMPA0820, add the following 
16 cards : 


IF(DELTR.GE. DDELO) GOTO 149 
WR1TE(MWR1TE,2000) DELTR 
2000 FORMAT (/,36H DELTR LESS THAN 
2 ****** * ,/) 


DO 201 0 1 = 1 , IK 
DO 201 0 d=1 ,NF 
IF(IFLAG(I,J).NE.2) 
I BIG= I 
dF=d 

PAL=RPCC(I,d) 

TMIN = 0 .0 
I FLAG ( I , d )* 1.0 
GOTO 1 00 
2010 CONTINUE 
GOTO 1 50 
149 CONTINUE 


DDELD AND EQUAL TO 
GOTO 2010 


IMPA0821 

IMPA0822 

.D15.6.13H *** IMPA0823 

IMPA0824 
IMPA0825 
IMPA0826 
IMPA0827 
IMPA0828 
IMPA0829 
IMFACS2A 
IMPA082B 
IMPA082C 
IMPA082D 
IMPA082E 
IMPA082F 
IMPA082G 


In Subroutine IMPACT, modify statement IMPA0950 to read: 


2TMIN, LNTMIN, RPC , dF , IFLAG ,H , FH , NF, AY , AZ , RPCC) 


IMPA0950 


In Subroutine IMPACT after statement IMPA0990, insert the following 
5 cards: 


C 

C 


IF(LNTMIN.EQ.O) CALL UPDATE( 1 . 0D0 , TU , TW, VY , VZ , TFCGU ,TFCGW, TALFA. IMPA0991 
* VELFU , VELFW , VELFA , DELTR , IKK,N F, ICP , AY , AZ) IMPA0992 

... IMPA0993 

IF NO ACCEPTABLE CONTACT TIMES HAVE BEEN FOUND, GO TO END OF ROUTINE IMPA0994 
IF( LNTMIN. £Q .0)GO TO 150 IMPA0995 
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In Subroutine IMPACT, remove statements IMPA1030, IMPA1040, and 
IMPA1050. 

In Subroutine IMPACT after statement IMPA1080 add: 


JJJ-1 


IBPA1 085 


In Subroutine IMPACT after statement IMPA1330, add these 13 cards: 


IFJNPP. GT.O .AND. JJJ.EQ. 0) GO TO 4C00 
GO TO 4 02 0 

4000 WRITE (M WRITS, 40 10) PHAX 

401") FORMAT I/, 10 £1 H *) .‘WARNING! PENETRATION HAS OCCORED 
•CORRECTION. ‘/1U ,10 (1H*),* OVERRIDE ENACTED. PHAX = 
TT=TB+DELTE 
TMIN=DEETR 
«JF=J FI 
P At=P AL 1 
IBIG-IBIG1 
IFLAG{IBIG, JF)= 1 
GO TO 100 

4020 CONTINUE 


IMPA1331 
IMPA 1332 
. 1MPA1333 
WITHOUT IMPACT IMPA 133 4 
*, 315.6) IMP A 1335 

IK PA 133 6 
IMPA1337 
IMPA 1338 
IMPA1339 
IMPA 133A 
IMPA133B 
IMPA 133C 
IMPA133D 


In Subroutine TCONT, modify statement TCON0020 to read: 
2L NT HI N.RPC , NFTMIN.I FL AG ,H, FH, N F, AY,AZ,RPCC) 

In Subroutine TCONT, modify statement TCON0070 to read: 
2YELFW {1 ) ,IFLAG (5 1 ,6 ) ,H ( 1) ,FH ( 1) ,RPCC (51,6) 

In Subroutine TCONT, delete statement TCON1090. 

In Subroutine TCONT after statement TCON1120, add: 

IF (TM.GT.XI) IF1AG (INUM ,IF) = 0 

In Subroutine TCONT after statement TCON1310, add: 

IF (P2CT »LT. Q.0.OR.P2CT.GT.P 1) IPLAG {ENUfl.IF) *0 


TCQN0020 


T CON0O70 


TOOK 1 125 


XCON1315 
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In Subroutine TCONT after statement TCON1320, add: 

RPCC (LNUa,IP) rn P2CT/P1 TCON1325 

In Subroutine TCONT after statement TCON1360, add: 

IP (TH.GT.TJJIR) GOTO 100 

TOOK 1365 


The message 

******WARNING! PENETRATION HAS OCCURRED WITHOUT IMPACT CORRECTION. 

***** * 0 VERR IDE ENACTED. PMAX = some number 

may printout just before the line which specified the IMPACT NO., TIME, CYCLE, 
ELEMENT, FRAG, and distance for an impact. This message indicates that an 
impact occurred numerically so close to the end of a timestep that the exact 
time-of- contact solution failed to detect the impact. However, a geometric 
penetration was detected and the impact interaction calculations were performed 
at the end of the time cycle. 

After the line which describes the impact is printed, a message which 

says 

AINT = a negative number NO IMPACT — LEAVING IMPCTE 

may print. This message indicates that a pseudo impact was erroneously 
detected by Subroutine IMPCTE and should be ignored. AINT is the relative 
velocity between the fragment and the ring, a negative value of aint 
indicates separation. This situation may occur at the initial point of 

fragment rebound. When this message prints, no impact interaction calculations 
have been performed. 
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SECTION 4 


MODIFICATIONS OF THE INPUT AND OUTPUT 
INSTRUCTIONS PERTINENT TO CIVM-JET 4C 


The input and output instructions given in Ref. 1 apply also to 
CIVM-JET 4C except for the following modifications (numbered 1 through 7) : 

1. Cards 7AA for the static "stress-strain" pairs for the mechanical sub- 
layer material model for the "main structure" (called structure 1) are 

EPS (1,1) and SIG(1,1) for the first coordinate pair 

EPS (2,1) and SIG(2,1) for the second coordinate pair 


are now defined by coordinate pairs denoting piecewise linear 
fits to: 


T = SIG versus £* = EPS 
u u 

o 

where = o^,(l+E ) = Kirchhoff uniaxial stress (subscript "o" 

o denotes static conditions) 

£* = £n(l+E^) = uniaxial logarithmic strain 

p 

and U E = — = engineering stress of a uniaxial static tensile 
o 

test specimen? P is the applied load and A is 

o 

the pre-test cross-sectional area of the specimen 

E^ = the measured axial relative elongation of the 
uniaxial test specimen 

= change in gage length 
original gage length 


= output which strain gages or extensometers 
can provide 
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In preparing the uniaxial test data in T u vs. e* form, the data in 
the strain region where necking occurs (that is, beyond the peak 
in a E = P/A q ) should be modified appropriately to "correct for 
necking". Various schemes for making such corrections have been 
developed. See, for example, the procedure and correction factor 
proposed by Bridgman [4] based upon extensive experimental and 
theoretical work. For more information on necking, see the book by 
Lubahn and Felgar [5] . Recent work on computer simulations of 
tension tests of ductile metals is reported by Norris et al. [6] and 

by Saje [7] . An excellent recent survey article on this subject was 
prepared by Hutchinson [8] . 

One approach to approximate the uniaxial behavior beyond the 
incipient necking condition (peak in 0 ) is to assume a straight 
line fit between that point and the rupture condition which can be 
characterized by the load at rupture and the cross-sectional 
area A f of the specimen measured at the rupture station after the 
test (ignoring any elastic recovery) . Hence, one can estimate the 
true stress at rupture as (c^) = P f /A f ; since at the associated very 

large strains, the material may be regarded as behaving in an 

incompressible fashion one may use the approximation p/p = 1. 

° 

Accordingly at rupture (failure) (x ) = (a ) since 

p ? ? H»f 

a T = A = P~ A~ * l4E u ) = D~ T u * The associated (£*) = Zn[l+(E ) ], 

O O *0 u -F u p 

P A A r 

where (E u ) is given by [1+(E ) ] = -2-° = . Finally, the 

f u f p A f A f 

"corrected" value to be used for (X ) is called (x ) and may be 

U *f ^fc 

computed, for example, by using Bridgman's [4] correction factor 
by (see Eq. 5-8 of Ref. 5) : 
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where 


-a = radius of the ( assumed to be circular) rupture cross 
section 

R = lateral final radius of curvature of the tensile test 
specimen at the rupture station. 

Bridgman [4] presents data plots (from extensive experiments) from 
which one can determine the ratio a/R from a knowledge of A q /A^. 
Other correction alternatives may be found in Refs. 5-8. 


2. The strain-displacement relation (called Type F for finite strain) 

from which the circumferential strain £ (or y of Ref. 1 hereafter 

2 11 
called y ) for the ring or beam structure is computed [2(p. 339) ; 

3 (Section 4)] : 




Ci +2 to 



(4.2) 


The through-the-thickness normal strain y 


3 

3 


is given by 



(4.3) 


since the material behavior is approximated as being incompressible, 

and it is assumed that the across-the-width strain of the ring is 
1 °2 

zero; that is, = 0. In the above, y^ represents the circumferen- 
tial (membrane) strain at the reference axis C = 0; £ refers to 
the location of particles before deformation (superscript "o"); 
and K denotes the "change of curvature": 

t* = X + jr7C + tV 2 (4 - 4) 
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(4.5) 



Figure 1 illustrates the geometry, coordinates, displacements, and 
generalized displacements for a curved-beam finite element, including 
the reference-surface displacements v and w as well as the quantities 
<J> and x* Note that indices 1, 2, 3 correspond to directions £, n, £. 
The present procedure used in CIVM-JET 4C accounts for finite strains 
and arbitrarily large rotations, and also accounts approximately for 
thickness changes accompanying finite strains in the plastic range, 
unlike CIVM-JET 4B which is valid only for small strains 3 and 
moderately small rotations . 


3. 


In CIVM-JET 4C, the uniaxial vs. e* data are represented by 
piecewise-linear segments which in turn are used in the mechanical- 
sublayer material model. Time-dependent plasticity is employed 
whereby it is assumed that strain rate affects the yield stress of 
the kth elastic, perfectly-plastic sublayer k x y according to: 



(4.6) 


where d and p are material "strain rate" constants, k r y is the 

• u o 

static yield stress of the kth sublayer, the rate of deformation is 


and 




(l +2 tfjf) 


(4.7) 


Ji 



All 

At 


(4.8) 


a: Defined by 1 + 2 y^ s 1, where is the circumferential strain, 
b: Defined by sinQ * 0 and cos0 s 1, where 0 is the angle of rotation. 
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(4.9) 


AlT, — y at time t 
2 Z 


minus y. 


at time 


(t - At) . 


The elastic and plastic energies as calculated by the ENERGY subrou- 
tine in CIVM-JET 4B and in CIVM-JET 4C are meaningful only for small 
strains. However, the "Total Strain Energy” calculated by the STRESS 
subroutine and printed in CIVM-JET 4C just after the Elastic and 
Plastic Energies is valid for finite (and small) strains. 

The treatment of the continuity between the "main structure material" 
and the "branch material" as presently handled by the BRAN subroutine 
in CIVM-JET 4B and in CIVM-JET 4C is incorrect for the case in which 
the elastic modulus of the branch material differs from that of the 
main- structure material (both for small and finite strains) . 

The behavior for the rotational springs of subroutine QREM (in both 
CIVM-JET 4B and CIVM-JET 4C) is valid only for the case of small 
strains and moderately small rotations at the locations of the 
rotational springs. 

To permit the user to label the printed output with a header line 
describing the case being run, one additional input card is required. 
Before Card #1 of the input, insert a Card #0 which should contain 
a description of the case in columns 1 through 80 . This card will 
be the first line of the printed output. If the user does not want 

a header line, this card should be left blank but it must be provided. 

In front of card 1, inserted the following card (called Card 0) : 

NEDI(J) for J = 1,2,3 ... 20 with format 20A4 


where 

the array NEDI is an 80 character heading. 


In the CIVM-JET 4C MAIN program, these 4 cards should be inserted 


before card 

MAIN0770: 


5555 

READ (MREAD, 5556) (NEDI(J) ,J=1,20) 

MAIN0761 

5556 

FORMAT (20A4) 

MAIN0762 


WRITE (MWRITE, 5557) (NEDI(J) ,J=1,20) 

MAIN0763 

5557 

FORMAT (1H1,20A4) 

MAIN0764 

and card MAIN0770 should be changed to read: 



READ (MREAD, l)B(l) ,DENS(1) ,EXANG, 

MAIN0770 


Finally, the reader is reminded that for both CIVM-JET 4B and CIVM-JET 4C, 
an assumed-displacement type of finite element is used. A cubic polynomial 
in the circumferential reference-surface coordinate r) is used for each of the 
reference-surface displacements v and w (see Fig. 1); this is called a cubic- 
cubic (CC) displacement element. Although this is a "compact" and satisfactory 
element, it accommodates at the element nodal junctions membrane strain 
continuity, but the bending contributions to the strain are not continuous. 
Higher order elements which provide continuity in (a) both of these strain 
contributions and (b) the n-direction derivative of each of these strain 
contributions have been investigated; the results are reported in Ref. 9. 
Although improved strain predictions result from using these higher-order 
elements, the storage and computing requirements also increase. Consequently, 
it was concluded [9] that it is more effective and efficient to employ a 
greater number of structural-model degrees of freedom (DOF) and CC elements 
than to use a comparable number of DOF's together with a smaller number of 
higher-order elements. Hence, the use of CC elements has been retained in 
both CIVM-JET 4B and CIVM-JET 4C. 
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SECTION 5 


ADDITIONAL USER INSTRUCTIONS FOR CIVM-JET 4C 

The CIVM-JET 4C computer program for the finite-strain analysis of 2-D 
transient structural responses of complete or partial rings and beams subjected 
to fragment impact has been stored on tape as a series of individual files. 

The attached Appendix A describes in detail which subroutines are found in 
these files. All references to the CIVM-JET 4C program are made assuming that 
the user has a copy of Ref. 1: NASA CR-134907 (ASRL TR 154-9) which serves 

as a user's guide to (1) the CIVM-JET 4B computer code and (2) the CIVM-JET 4C 
computer code with the use of the modified input instructions given in the 
present report. All of the subroutines for CIVM-JET 4C are the same as those 
for CIVM-JET 4B except for Subroutines MAIN, ELMPP, ENERGY, IMPACT, PRINT, and 
STRESS ; these latter subroutines have been given new sequence numbers . A 
complete FORTRAN IV source listing of the CIVM-JET 4C subroutines accompanies 
the tape which is available from MIT under a copyright licensing agreement. + 

This is usually a 9- track tape with the program stored in card- image 
format: IBM EBCDIC code, 80 characters to a record, 1600 characters to a block, 
and 1600 bits per inch. The tape does not contain any labels; therefore, a 
no- label format must be used when accessing the tape. A magnetic tape with 
other characteristics such as 7-track, ASCII code, 800 bits per inch, etc. 
could be provided; the desired format should be requested. 

This program was formulated and run on the IBM 370/168 computer facility 
at MIT. The tape is set up to be compiled and run immediately on an IBM 
machine by either direct access or by being transferred to disc storage. 
However, if a non-IBM system is used, care should be taken to make sure that 
the program is compatible with the operating system. The CIVM-JET 4C program 
assumes that the operating system uses code 5 for a reader, code 6 for a line 
printer, and a code 7 for a card punch. If your computer uses a different 
convention, it will be necessary to obtain a card deck from the tape and 

+ To obtain CIVM-JET 4C, please contact Prof. E.A. Witmer, Room 41-219, 

MIT, Cambridge, Mass. 02139. 
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change the code numbers in the MAIN program; see page 62 of Ref. 1 for the 
cards that would need to be changed. Note also that the symbol @ has been 
employed to denote a continuation card; if this is unacceptable for your 
computer facility, change this appropriately. 

To check the operation of this tape, it is suggested that the user run 
the two examples (6.1 and 6.2) listed in Ref. 1. input data for each of the 
two examples have been included on the tape for this purpose; outputs from 
CIVM-JET 4C for both of these examples are supplied to acquirers of the tape/ 
Because of the large number of computational steps used in the solution of 
each of these problems and because of the ever-changing IBM systems and 
subroutine libraries, the user should be forewarned that exact duplication 
of these example runs with the results attached hereto probably will not be 
obtained. Instead, slight variances will be detected because of roundoff; 

Table 2 illustrates a sample of the differences in certain solution quantities 
for each of problems 6.1 and 6.2 between IBM 370/168 calculations at MIT and 
UNIVAC calculations carried out at NASA-Lewis [10] . The purpose of this 
exercise is to be certain that all of the options of the CIVM-JET 4C computer 
code are compatible with the user's computer system. 

In order to run the program, any one of three methods may be utilized: 

(1) Direct Access: From Appendix A find the files that are needed for 
a particular run and compile each one, concatenating the compiled 
codes into one block in your computer code. 

(2) Using Disc: The same as above, except after compiling the files, 

they are stored on disc in object format form for future use without 
further compilation. Concatenation will still be necessary. 

(3) Card Format: Having the programs punched onto cards directly from 

the tape, and use the cards to run the program. This is necessary 

if any changes are to be made in any of the subroutines. The subrou- 
tines that are to be changed must then be punched onto cards, but the 

rest of the program can still be accessed by using methods (1) or (2) 
listed above. 


See footnote on page 22. 
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TABLE 1 


FEATURES AND CAPABILITIES OF THE CURRENT CIVM-JET 4B AND 
CIVM-JET 4C COMPUTER PROGRAMS 


Feature 

CIVM-JET 4B 

CIVM-JET 4C 

Spatial Approximation 
Finite Element 

X 

X 

Timewise Central 

Difference Operator 

X 

X 

Ring Geometry 

Complete Ring 

X 

X 

Partial Ring 

X 

X 

Initial Configuration 

Circular 

X 

X 

Arb. Curved 

X 

X 

Constant Thickness 

X 

X 

Variable Thickness 

X 

X 

Single Layer 

X 

X 

Boundary Conditions 

Ideally Clamped 

X 

X 

Hinged Fixed 

X 

X 

Free 

X 

X 

Clamped-S liding 

X 

X 

Other Support Conditions 

Distributed Elastic Foundation 

X 

X 

Point Elastic Springs 

X 

X 

Structural Branch 

X 

X 




TABLE 1 (CONTINUED) 



CIVM-JET 4B 

CIVM-JET 4C 

Material 



Single Material 



for Each Branch 

X 

X 

Homogeneous 

X 

X 

Initially Isotropic 

X 

X 

Temperature Independent 

X 

X 

EL 

X 

X 

EL-PP 

X 

X 

EL-LSH 

X 

X 

EL-SH 

X 

X 

EL-SH-SR 

X 

X 

Small Strain 

X 

X 

Finite Strain 

- 

X 

Stimulii 



Impacting Fragments 



Single 

X 

X 

Multiple 

X 

X 

Friction 

X 

X 

Elastic (e=l) 

X 

X 

Inelastic (e=0) 

X 

X 

Intermediate 0<e<l 

X 

X 

Deflections : Bernoulli-Euler 



Type Only 



Small 

X 

X 

Large 

X 

X 

Strains 



Membrane 



Small 

X 

X 

Finite 

X 

X 

Bending 



Small 

X 

X 

Finite 

- 

X 
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TABLE 1 (CONTINUED) 


Strains (Continued) 
Overall 
Small 
Finite 

CIVM-JET 4B 
X 

CIVM-JET 4C 

X 

X 

OUTPUT INFORMATION 



At Selected Times 



Energy/Work Type and Amount 



Structure Kinetic Energy 

X 

X 

Structure Strain Energy 

X 

X 

Elastic: Small Strain 

X 

X 

Plastic: Small strain 

X 

X 

Strain Energy of Elastic Restraints 

X 

X 

Total 



Small Strain 

X 

X 

Finite Strain 

- 

X 

Fragment Kinetic Energy (Each) 



Translational 

X 

X 

Rotational 

X 

X 

Total 

X 

X 

Circumferential Strains 



Inner Surface 

X 

X 

Midsurface 

- 

X 

Outer Surface 

X 

X 

At Gaussian Stations 

X 

X 

At Nodal Stations 

X 

X 

At Additional Stations 

X 

X 

Location where Prescribed 



Strain is Exceeded 

X 

X 

Support Reactions 

X 

X 
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TABLE 1 (CONCLUDED) 


Time, Location, and Value of 
Largest Strain Reached During Run 

CIVM-JET 4B 

X 

CIVM-JET 4C 

X 

CAPACITY INFORMATION 



Max. No. of Finite Elements* 

50 

50 

Max. No. of Fragments* 

6 

6 


These limits can be circumvented by altering the dimensions of appropriate 
program variables (see each source reference) . 
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TABLE 2 


COMPARISON OF CIVM-JET 4C RESULTS FROM USING DIFFERENT COMPUTERS 


System 

Compiler 

Site 

Example" 1 " 

M 

STR 

Largest Added 
Pt. Strain 

Largest Nodal 
Strain 



Substructure 



i 

2 

i 

2 

i 

2 

IBM 370 

H 

MIT 

6.1 

.0596839 

.0129698 

.0179878 

.00184006 

.0648378 

.0142004 

UNIVAC 

ASCII 

NASA-Lewis 

6.1 

.0596839 

.0129094 

.0179878 

.00184006 

.0648371 

.0140876 

IBM 370 

H 

MIT 

6.2 

.154770 


.0968485 


.185669 


UNIVAC 

ASCII 

NASA-Lewis 

6.2 

.154863 


.0958849 


.185654 



These are the examples defined in Ref. 1 but with input appropriate for CIVM-JET 4C 
as noted in Section 5 of the present report. 




REFERENCE CONFIGURATION 



-15° < <J>. _ - <p. < 0 
-180° < d>. < 180° 

r i — 

<!> (n) = b Q + b x n + b 2 u 2 

R(n) = -0^/9n) -1 N0DE i; n = o 

h(n) = h. (i - ^-) + h node i+i: n = n. 

1 Tj. 1+J. T| . -L 


LOCAL SYSTEM CARTESIAN REFERENCE 


£» Hr C 

COORDINATES . 

Y,Z . 

- - GLOBAL 


• 


COORDINATES 

v,w ,if»,x 

DISPLACEMENTS 

Y/2 

- LOCAL 

A 

H 

A 

to 

• 

• 

• 

A 

00 

1 

ELEMENT GENERALIZED 
DISPLACEMENTS 

- 

COORDINATES 

q l *2 «3 q 4 _ 

v. w - X- 

X X r x A x 



v 3n r 

y = *1 + ” 

x an r 



FIG. 1 NOMENCLATURE FOR GEOMETRY, 

COORDINATES, 

AND DISPLACEMENTS OF 


A CURVED-BEAM FINITE ELEMENT 
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APPENDIX A 



TAPE 

FILES FOR 

CIVM-JET 4C 



The CIVM-JET 

4C program 

consists of 

the following 

main program and 23 

subroutines ; the 

file numbers 

associated 

with each subroutine 

are indicated 

File 

Subroutine 

File 

Subroutine File 

Subroutine 

1 

MAIN 

9 

FICOL 

17 

CUBIC 

2 

ASSEF 

10 

IDENT 

18 

IMPACT 

3 

ASSEM 

11 

MINV 

19 

IMPCTE 

4 

BRAN 

12 

OMULT 

20 

PENTRN 

5 

DINIT 

13 

PRINT 

21 

ROOT 4 

6 

ELMPP 

14 

QREM 

22 

ROTAT 

7 

ENERGY 

15 

STRESS 

23 

TCONT 

8 

ERC 

16 

TSTEP 

24 

UPDATE 


File 25 Input Data for Example 6.1 of Ref. 1; applicable to both 
CIVM-JET 4B and CIVM-JET 4C. 


File 26 Input Data for Example 6.2 of Ref. 1; applicable to both 
CIVM-JET 4B and CIVM-JET 4C. 

The program is written in double precision arithmetic. The number of 
memory locations required on the IBM 370/168 computer at MIT is approximately 
440,000 bytes; this includes locations for the MIT computer library subroutines. 


Subroutines 

Related to 

the Finite Element Solution Procedure: 


ASSEF 


ENERGY 

OMULT 


ASSEM 


ERC 

PRINT 


BRAN 


FICOL 

QREM 


DINIT 


IDENT 

STRESS 


ELMPP 


MINV 

TSTEP 

Subroutines 

Related to 

Fragment-Ring- 

-Impact Interaction: 



CUBIC 


PENTRN 

TCONT 


IMPACT 


ROOT 4 

UPDATE 


IMPCTE 

31 

ROTAT 







